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In the statistical mechanics approach introduced by 
S.F. Edwards to describe granular media in their me- 
chanically stable states, the granular entropy is defined 
as the logarithm of the number of mechanically stable 
configurations N grains can be arranged in a volume V. 
In order to calculate this entropy, it is necessary to intro- 
duce a volume function W expressing the total volume 
as a function of the position and coordination of the N 
grains. Recently Blumenfeld and Edwards [l| have pro- 
posed such a volume function in two dimensions. This 
volume function is constructed by first partitioning the 
system into voids defined from the grain contacts, and 
then by partitioning the voids into quadrilaters defined 
from the grains centers and the position of the contact 
points. The volume of the system is equal to the sum of 
the volumes of the voids, i.e. to the sum of the areas of 
the quadrilaterals. Here I show that the proposed parti- 
tioning of the voids into quadrilaterals is not exact, and 
therefore that the introduced volume function is not a 
volume function as claimed. 

The partition of a void into quadrilaterals works when 
the void is convex, as that shown in Fig. 1 of Ref. [![. 
Some non-convex voids, which has been overlooked in 
Ref. are instead not well partitioned. In order to 
show that the presence of these non-convex voids is rele- 
vant, I have investigated the voids of mechanically stable 
(under gravity) configurations of 2D polydisperse gran- 
ular systems, generated by simulating a sedimentation 
process via standard soft-core Molecular Dynamics sim- 
ulations 0. Fig. Q] shows (the central region) of one of 
the investigated configurations, together with the parti- 
tion of the system into voids. In this configuration about 
20% of the voids are non convex, and 5% of the voids (the 
shaded ones) are not well tesselleted by the quadrilaterals 
introduced in Ref. 

Two parameters are of interest to evaluate the effect of 
non-convex voids on the analysis performed by Blumen- 
fcd and Edwards. One is the fraction / of voids which 
are not well tesselleted, the other is the relative error in 
the estimation of the volume, rj = (W — V)/V, where V 
is the volume occupied by the grains, and W its estima- 
tion. Figure 2a and 2b shows the dependence of / and 
r\ on the polydispersity of the system [4j, for /z = 0.2, 0.4 
and 0.6. The magnitude of r) is understood by consid- 
ering as an example a monodisperse system of volume 
V and area fraction (f>, with total volume V — nv s /<p, 
n number of particles and v s area of a particle. In this 
case (fr r ip ~ 0.84 < 4> < 4>crystai — 0.90, and a very crude 
extimation of the total volume as W = nv s /(f> rn , with 
<t>m = {4>ri P + 4> 'crystal) /2, gives a relative error rj < 4%. 




FIG. 1: The central part of a mechanically stable (under 
gravity) configuration of a polydisperse system is tessellated 
in voids according to the procedure of Ref. flj. The shaded 
non-convex voids are not well tessellated by the quadrilaters 
defined in Ref. [J]. 
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FIG. 2: (a) The fraction of voids / not well partitioned by 
the quadrilaterals introduced in Ref. [l{; (b) The relative error 
made when measuring the volume with the volume function 
introduced in Ref. 



These results indicates that / and rj increase both with 
p and with /i, and that for spherical grains the volume 
function introduced by Blumenfeld and Edwards works 
only at first crude approximation, even in the absence 
of polydispersity. The presence of non-convex voids will 
be certainly more relevant when non-spherical grains are 
considered. 
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